Measurement of the velocity of neutrinos from the CNGS beam with the Large 

Volume Detector 
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We report the measurement of the time-of-fiight of ~ 17 GeV on the CNGS baseline (732 km) 
with the Large Volume Detector (LVD) at the Gran Sasso Laboratory. The CERN-SPS accelerator 
has been operated from May 10th to May 24th 2012,with a tightly bunched-beam structure to allow 
the velocity of neutrinos to be accurately measured on an event-by-event basis. LVD has detected 
48 neutrino events, associated to the beam, with a high absolute time accuracy. These events allow 
to establish the following limit on the difference between the neutrino speed and the light velocity: 
-3.8 ■ 10~ 6 < (v v -c)/c < 3.1 ■ 10~ 6 (at 99% C.L.). This value is an order of magnitude lower than 
previous direct measurements. 

PACS numbers: 06.30.Gv,14.60.Lm,29.40.Mc 



Introduction.- Cosmological measurements [lj provide 
stringent limits on the sum of neutrino masses. Even 
assuming the heaviest neutrino eigenstate, the expected 
relative deviation from speed of light is lower than 10~ 20 
for 10 GeV neutrinos. Nevertheless, in the past, theo- 
ries allowing some or all neutrinos to have apparent ve- 
locities different than the speed of light, c, (see e.g. 0) 
have been proposed. A stringent limit at E ps 10 MeV, 
\v v — c|/c < 2 x 10~ 9 , has been obtained from the obser- 
vation of SN1987A electron anti-neutrinos Q. At higher 
energies (E > 30 GeV), the deviation has been tested 
down to L - c\/c < 4 x 10" 5 Q. The MINOS collab- 
oration [5J has performed a neutrino time-of-flight (tof) 
measurement on a w 735 km baseline and with a beam 
with average energy < E >=3 GeV. The MINOS result 
is quoting (v u - c)/c = (5.1 ± 2.9) x 10~ 5 at 68% C.L.. 
Recently the OPERA collaboration has reported a tof 
measurement for muon neutrinos from the CNGS (CERN 
Neutrinos to Gran Sasso) beam. The reported evidence 
of a superluminal propagation of was subsequently at- 
tributed to a technical problem Q . This explanation has 
been also confirmed by the results of another experiment 
in the Gran Sasso Laboratory [§]. 

The Large Volume Detector (LVD), in the INFN Gran 



Sasso National Laboratory (LNGS), at the average depth 
of 3600 m w. e., is a 1 kt liquid scintillator detector whose 
major purpose is monitoring the Galaxy to study neu- 
trino bursts from gravitational stellar collapses [9j. It 
has started operation in 1992, and since 2006 it has been 
acting as a far- monitor of the CNGS beam 1(J 11 1 . LVD 



is sensitive to neutrino interactions with protons and car- 
bon nuclei in the liquid scintillator and with the iron of 
the detector structure. Muons, produced by charged cur- 
rent interactions of muon neutrinos in the rock, can also 
be detected and are responsible for the bulk of CNGS 
events in LVD. We report a measurement of the neutrino 
velocity obtained through the detection of from the 
CNGS beam. We show that the neutrino speed is com- 
patible with c, its deviation being < 3.8 x 10~ 6 at 99% 
C.L.. The sensitivity of this measurement is by one or- 
der of magnitude better than that of previous ones in a 
similar energy range. 

The detector. - LVD consists of 840 scintillator coun- 
ters, 1.5 m 3 each. The array is divided in three identi- 
cal "towers" with independent high voltage power sup- 
ply, trigger, data acquisition and absolute clock (ESAT 
Slave RAD100) connected, through a 8 km optical link, 
to the Master clock (ESAT RAD100) located in the ex- 
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ternal buildings of the LNGS. Each tower consists of 35 
"modules" hosting a cluster of 8 counters. Each counter 
is viewed from the top by three 15 cm photomultiplicr 
tubes (PMTs) FEU49b or FEU125. 
LVD standard electronics is described in detail in 12 1 13j |. 
The trigger condition for each tower is the three-fold 
coincidence of the PMTs of any of its counters, corre- 
sponding, on average, to the energy threshold Eh — 4 
MeV. The energy released in every counter is measured, 
through 12 bit ADCs, with a mean resolution cfe/E ~ 
15% at 10 MeV. The arrival time is measured with a 12.5 
ns granularity. A high stability Citrine crystal (40 MHz) 
oscillator supplies the general clock for the whole exper- 
iment. One millisecond after the occurrence of a trigger, 
the memory buffers, containing the energy and time in- 
formation of all detected signals, are read out together 
with the associated time stamp from every tower's abso- 
lute clock (ES AT Slave) . The absolute time accuracy of 
LVD is of the order of one microsecond. This is better 
than what needed to search for coincidences among differ- 
ent neutrino telescopes generated by a gravitational col- 
lapse and to guide the search for possible gravitational- 
wave signals 14 1. 

In spite of its limited absolute time resolution, LVD has 
detected, at the end of 2011, 32 neutrino events from the 
CNGS bunched-beam test [l5j. The difference between 
the neutrino time-of-flight toj v and the expected one at 
speed of light tof c has been found to be: 



5t=tof v -tof c = 3.1 ±5.3 stot ±: 
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To perform this measurement, we have carefully re- 
calibrated all LVD components and time delays in order 
to improve time accuracy. The uncertainties are due to 
the absolute time resolution; to the fluctuations in the 
time response of the different detectors (average values 
of transit times were used when not known) ; and to the 
uncertainty in determining the event position. This pre- 
liminary measurement has shown that, within the uncer- 
tainties, the speed of CNGS neutrinos is compatible with 
that of the light. Moreover, it has helped identifying the 
detector characteristics to be improved for a more accu- 
rate measurement of the neutrino velocity. 
Setup upgrades. - To allow a very accurate neutrino 
tof measurement, from May 10 th to May 24 th , 2012, 
the CERN-SPS accelerator has been operated with a 
new beam structure. This structure was made of four 
batches separated by about 300 ns, with 16 bunches per 
batch, with a narrow width of ~ 3 ns each, separated 
by 100 ns. The batched structure is known through the 
waveform obtained from the Beam Current Transformer 
(BCT) (see fig. [2j grey lines). 

For this occasion, a new High Precision Time Facility 
(HPTF) has been designed by the Borexino collabora- 
tion [16[ and installed in the external buildings of the 
LNGS. A Scptcntrio PolaRx4 GNSS receiver, synchro- 
nized with the 10 MHz frequency of a GPS disciplinatcd 



Rubidium clock, provides GPS time and XPPS (10 Hz) 
output signal. The HPTF is equipped with high pre- 
cision (50 ps) Time Interval Counters (TIC) Pendulum 
CNT-91, to which the triggers of the different LNGS ex- 
periments can be connected. 

In view of the new neutrino beam, we have modified a 
subset of the LVD counters, to improve their timing per- 
formances. We have chosen 58 of them (Super-Set, see 
fig. [T]) to maximize the acceptance with respect to CNGS 
neutrinos while minimizing the number of detectors to be 
modified. From data taken since 2006 from the CNGS 
beam (and in agreement with Monte Carlo simulations 
[Tol ]) we have measured that the Super-Set counters are 
involved in ~ 40% of the CNGS events detected by LVD, 
while representing only less than 7% of the whole array. 
On one hand, to avoid time fluctuations in the trigger for- 
mation at the single-counter level, we have modified the 
cabling of the PMTs, by delaying only one of them. The 
change guarantees that the 3-fold coincidence among the 
PMTs in every counter is always formed due to the same 
tube. On the other hand, to perform a measurement 
of the transit time in each counter, we have equipped 
them with a LED system. The transit time, denominated 
5 lvd, is the time between the light generation inside a 
counter and the formation of the trigger. By means of 
the LED system, we have measured the behavior of the 
transit time versus energy for each counter. By varying 
the LED intensity we have simulated different energy re- 
leases in each of them, and measured the transit time 
versus energy. For each counter, we have best fitted the 
results with a power-law of the form: 



5 lvd (E m ) — Pi ■ E„ 



(2) 



where E m is the energy release measured by the ADC, 
while PI, P2 and k are free parameters of the fit, for each 
counter. 

The trigger of the array has been upgraded too. An 
independent fast trigger logic has been implemented by 
extracting the coincidence of the triggered counters from 
discriminators, with the trigger sent to the HPTF. The 
trigger has been connected to one of the TIC, this being 
stopped by the XPPS signal. This provides a high preci- 
sion time difference, Stic, between the LVD trigger and 
the absolute GPS time. Thanks to this improvement, 
the absolute time accuracy of LVD is of the order of few 
nanoseconds. 

Finally a new, indipendent, high precision geodesy mea- 
surement has been performed. The CNGS-LVD base- 
line, namely the distance between the center of the 
BFCTI.400344 Beam Current Transformer (BCT) de- 
tector at CERN and the LVD Super-Set upstream en- 
try wall (taken as the LVD reference), is found to be 
731291.87±0.04 The corresponding time-of-flight 

at speed of light, when including the 2.2 ns contribution 
due to Earth's rotation, is tof c = 2439329.32 ± 0.13 ns. 
The Measurement. - The time-of-flight of neutrinos, fo/„, 
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FIG. 1: Schema of LVD counters. The darkened (blue) ones 
represent the Super-Set, see text. 



is the difference between the absolute time at which the 
event triggers a counter in LVD (t u ) and the absolute 
time at which the proton bunch crosses the BCT inten- 
sity monitor at CERN (t p ), both expressed in the LNGS 
reference time, i.e., tof v = t v —t p . 
t v is given by: 



t 
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(3) 



'•xpps i s the absolute GPS time stamp of the first XPPS 
signal after the trigger. 8h is a correction for the path 
difference due to neutrino events hitting counters at 
different heights (the angle of incidence of the CNGS 
beam with respect to the horizontal plane is « 3 deg). 
^■const — 180.8 ± 1.5 ns sums up all constant delays in 
LVD. 5 lvd and <5tic have been previously introduced. 
They account for the transit time correction of each 
counter and the delay between the trigger and the ab- 
solute GPS time, respectively. 
t p is given by: 
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where tncA is the time stamp of the kicker proton extrac- 
tion signal (start of the BCT waveform digitalization); 
5 w f is the delay of the first proton bunch with respect 
to the beginning of the waveform acquisition (tncA)', 
AgVRN = 9521.1 ± 2.0 ns sums up all the constant de- 
lays at CERN and SlnGS 1S the time difference between 
the LNGS and CERN reference systems \v^. 
LVD has been fully operational during the May 2012 
bunched-beam run. For each detected event in that pe- 
riod, we have determined the difference between the time 
of detection in LVD and that of the kicker extraction sig- 
nal at CERN, after accounting for the CERN-LVD base- 
line, i.e., dT batch = (U - t p ) - tof c . 5T batc h does not 



include yet the identification of the neutrino bunch orig- 
inating the event. To check the detector operation, we 
have looked for coincidences in the whole LVD in a rather 
wide time window, namely \STf, a tch\ < 100 [is. We have 
found 190 events in total, consistent with the 1.89 • 10 17 
protons on target (p.o.t.) delivered during the bunched- 
beam run [l(| ■ 

For the measurement of the neutrino velocity, we have 
used only events involving at least one of the Super-Set 
counters. We have found 79 of them out of 190, i.e., 40% 
as expected (see previous section). Four of them are ac- 
tually not usable for an accurate time-of-flight measure- 
ment: one has no complete information in the CERN 
database, while for three of them the ADC was malfunc- 
tioning. To limit the sources of systematic uncertain- 
ties, we have applied quality cuts to the remaining 75 
events. First, because the number of photoelectrons for 
energy releases E<10 MeV is too low to guarantee that 
the counter is triggered by direct (i.e., fast) light, we 
have selected only events where E>10 MeV (ten excluded 
events). Second, we have required the non-saturation of 
ADC of the triggering counter, as the saturation would 
not allow a precise measurement of the Slvd term in 
(2) (seven excluded events) . Finally on the remaining 58 
events, we have tested different energy cuts, between 10 
and 100 MeV. While the mean value of the time-of-flight 
distribution remains constant for any cut, the r.m.s. de- 
creases as the threshold increases. It becomes stable for 
E>50 MeV, this value being our final choice for the en- 
ergy cut. The described cuts reduce the sample to 48 
events, that will be used for the final analysis. We note 
that among all the excluded events, there are five whose 
time-of-fiight significantly deviates from the mean of the 
distribution of the original 79 events. Those are five of 
the events for which E<10 MeV. Their number is con- 
sistent with that expected from background. Indeed, the 
background frequency of the Super-Set, mainly due to 
the rock radioactivity, is f^h = 1.1 s _1 . With n ~ 25000 
extractions, the expected number of background events 
is N = nf bk 2 • 10~ 4 w 5.5. 

The 48 selected events have been individually inspected: 
almost all of them show the presence of a muon track, 
as expected from a Vp charged current interaction. We 
show the distribution of the measured ST^atch in figure 
[2] (black lines) where the sum of the associated digitized 
waveforms of the beam intensity monitor are also shown 
(grey lines). The beam structure is clearly identifiable, 
and every LVD event can be associated to the closest 
beam-waveform peak. For each event we then calculate 
the time difference, St, with respect to the maximum in- 
tensity of the peak. The distribution of St for the 48 
events is shown in figure [3] (black histogram) compared 
with the superposition of all the peaks of the waveforms 
correlated to detected events (grey curve). The positive 
tail in the grey histogram is an artifact of the transfer 
function of the BCT system [l^. This effect does not in- 
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fluence our measurement since time calibrations are per- 
formed with respect to the position of the maximum. The 
mean value of the measured distribution is: 

St = tof v - tof c = 0.9 ± 0.6 sta4 ns (5) 

The long-term stability of the experimental setup has 
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FIG. 2: Comparison of the ST values of the 48 selected events 
(black lines) with the summed waveforms of proton extraction 
(grey lines). The origin of time for the waveforms is given by 
the maximum of the first bunch. 
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FIG. 3: Distribution of St, the difference between the time- 
of- flight of neutrino and the time-of- flight at speed of light, 
for the 48 selected events (black histogram), compared with 
the superposition of the peaks of the waveforms correlated to 
detected events (grey histogram). 

been monitored during the measurement by triggering 
the LED in one of the counters of the Super-Set every 10 
minutes and measuring the delay of the entire chain (see 
figure |4} . The system, during the entire period, fluctu- 
ated with a r.m.s. =1.5 ns. The asymmetry of the distri- 
bution is due to PMT pre-pulses [lj| . 
The systematic uncertainties associated to the measure- 
ment of St are summarized in Table HI The baseline be- 
tween LVD and CERN is known with an uncertainty 
lower than 0.3 ns. In fact, neutrino-induced muons re- 
sponsible for events in the Super-Set can be generated as 
far as several hundreds of meters from LVD. This makes 
shorter the actual neutrinos baseline, while part of the 
distance is travelled by muons of different energy. This 
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FIG. 4: System stability during the measurement: every 10 
minutes the transit time of one Super-Set counter was mea- 
sured. 



Systematic uncertainties 


ns 


Error dist. 


Baseline (< 10 cm) 


< 0.3 


Gauss. 


Const, corrections at CERN 


2.0 


Gauss. 


BCT calibration 


1.0 


Gauss. 


Time link calibration 


1.1 


Gauss. 


GPS syncronization accuracy 


1.0 


Gauss. 


PMT transit time correction 


0.4 


Gauss. 


Absolute transit time calib. 


1.5 


Flat 


Optical fiber length 


0.5 


Gauss. 


Fluorescence lifetime 


0.6 


Exp. 


Total systematic 


± 3.2 





TABLE I: Sources of systematic uncertainty in the measure- 
ment. 



effect has been investigated through Monte Carlo and 
has a negligible impact on the measurement. Constant 
delays at CERN give an additional uncertainty of 2.0 ns. 
The uncertainty due to the system time inter-calibration 
[l|| has been measured and it is 1.1 ns. Also, the delay 
between the proton extraction time and the recording of 
the BCT waveforms by a digitizer has been measured, 
with an uncertainty of 1.0 ns. The variable correction 
Slvd (energy and counter dependent) introduces a fur- 
ther systematic uncertainty, 0.4 ns, evaluated through 
the propagation of the one associated to the energy mea- 
surement. Finally, the last three terms in Table U are 
related to the constant corrections in the absolute time 
calibration of the Super-Set counters. They account for 
the unknown time delay between trigger and LED light 
generation, 2.6±1.5 ns; the uncertainty associated to the 
optical fiber length, 0.5 ns, and the difference between 
the detector response to LED light and to particle ion- 
ization, 1.1 ±0.6 ns (the scintillator decay time is r=3.32 
ns By quadratically summing up all these contri- 

butions, we obtain a total systematic uncertainty of 3.2 
ns. 

Conclusions.- We have presented the measurement of the 
neutrino velocity with the LVD experiment, through the 
detection of muon-ncutrinos from the CNGS beam. Dur- 
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ing 10 days of bunched-beam, LVD has detected 190 
events in coincidence with the beam neutrinos. Among 
them, 79 have involved at least one counter in the Super- 
Set, an ensemble of 58 counters upgraded for precise tim- 
ing measurements. To limit the sources of systematic un- 
certainties, we have applied quality cuts to these events. 
The resulting 48 have been used to determine the time-of- 
flight of with < E > = 17 GcV on the CNGS baseline. 
The deviation from the time expected from propagation 
at the speed of light has been found to be: 

St = 0.9 ± 0.6 stat ± 3.2 sys ns (6) 

The corresponding 99% confidence limit on the speed of 
neutrino is: 

-3.8xlO~ 6 <K-c)/c<3.1xlO" 6 (7) 

These values are an order of magnitude lower than pre- 
vious direct measurements. 

By using the average neutrino beam energy, < E v > = 17 
GeV, a limit is found on the relativistic mass of neutrino: 

m„„ < 47 MeV/c 2 (99% C.L.) (8) 

Acknowledgements.- This measurement was made possi- 
ble by the effort of G. Di Carlo, S. Parlati and P. Spinnato 
from the LNGS Computer Center; G. Korga, from the 
Borexino Collaboration; P. Alvarez and J. Serrano from 
the CERN staff; G. Cerretto, V. Pettiti and C. Plan- 
tard from INRIM Torino; Hector Esteban from the Real 
Instituto y Observatorio de la Armada (ROA), San Fer- 
nando, Spain and the INFN Torino Electronics Division. 
We thank F. Vissani for discussions of great value to our 
work. 

Note added in proof. - Recently we received a draft con- 
cerning a new measurement of the A~^^ , i.e., the sum 
of constant delays at CERN [2l|. The new values is: 
A^f t N = 9522.4 ± 2.0 ns. By using it our result be- 
comes: 

St = -0.3 ± 0.6 stat ± 3.2 sys ns 

The difference with respect to the previous result is well 
within our given uncertainties. The corresponding limit 
on the speed of neutrino, at 99% C.L., becomes more 
stringent: 



and for the neutrino mass: 



m Vfl < 44 MeV/c 2 (99% C.L.) 
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